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We critically discuss recent suggestion to use long-range modes of charge (electric or baryon) 
fluctuations as a signal for the presence of Quark-Gluon Plasma at the early stages of a heavy ion 
collision. We evaluate the rate of diffusion in rapidity for different secondaries, and argue that for 
conditions of the SPS experiments, it is strong enough to relax the magnitude of those fluctuations 
almost to its equilibrium values, given by hadronic "resonance gas". We evaluate the detector 
acceptance needed to measure such "primordial" long-range fluctuations at RHIC conditions. We 
conclude with an application of the charge fluctuation analysis to the search for the QCD critical 
point. 



I. INTRODUCTION 

During the last few years significant attention has been 
attracted to the issue of event-by-event fluctuations in 
high energy heavy ion collisions. The original goal of 
such studies formulated in has been focused on equi- 
librium thermodynamical fluctuations at freeze- out. It 
was argued that long and intense final state interaction 
of secondaries makes heavy ion collision very different 
from hadronic collisions, in which dynamical fluctuations 
of quantum origin produce quite different "intermittent" 
behavior. Experimental data, such as obtained by NA49 
collaboration ||, have indeed revealed only very small 
Gaussian fluctuations of different quantities incompatible 
with earlier predictions of large non-equilibrium fluctu- 
ations, e.g. due to bubble formation during the phase 
transition. As discussed in [§-|6| in detail, extensive 
quantities like the total multiplicity, N, obtain contri- 
butions from final state interaction of secondaries due to 
resonances, and initial state fluctuations related to the 
number of participant/spectator nucleons. Furthermore, 
experimentally, for central collisions ((AN) 2 ) / (N) « 2, 
and both effects share about equally the responsibility 
for moving this number away from the Poissonian value 
((AN) 2 )/(N) = 1- As emphasized in ||fij, the intensive 
quantities, such as mean px, particles ratios etc. are de- 
termined entirely by the final state effects, and are well 
represented by the resonance gas. 

Other goals of the event-by-event fluctuation analysis 
have also been discussed. In particular, one of us us- 
ing the analogy between Big Bang phenomena and the 
"Little Bang" in heavy ion collisions, suggested to look 
for primordial "frozen QGP oscillations" , as it is done 
for anisotropic components of microwave background in 
the Universe. The idea was inspired by an observation 
that large-amplitude oscillations can be excited due to 



instabilities of parton counter-flows at the initial stage 

f . On the other hand, two simultaneous papers, by 
sakawa, Heinz and Muller [|l0| and Jeon and Koch |Tl| ] 
suggested a different set of primordial signals: equilib- 
rium charge fluctuations in QGP, which happen to be 
factor 2-3 smaller than in the hadron gas. 

The idea Jl(],[n| is based on the well known phe- 
nomenon of kinetic slowdown of fluctuation relaxation, 
provided sufficiently long-range harmonics of a conserved 
charge density are considered. If the relaxation time hap- 
pens to be shorter than the lifetime of hadronic stage of 
the collisions, the authors argue, the values of such fluc- 
tuations should deviate from their equilibrium (resonance 
gas) values towards their earlier, primordial values, typ- 
ical for QGP. 

Although this idea should work for parametrically 
long-range effects, whether a number of necessary con- 
ditions can indeed be fulfilled in realistic heavy ion col- 
lisions depends on the answers to the following ques- 
tions: (i) How fast is the relaxation due to final-state 
re-scattering in hadronic gas? (ii) How large the detec- 
tor acceptance should really be, in order to counter the 
relaxation? (iii) Is there a window of parameters, for- 
bidding resonance gas relaxation but allowing for it at 
the QGP stage? (iv) If not, what fluctuations should fol- 
low from non-equilibrium parton kinetics, at very early 
pre-QGP stages of the collisions? 

In this paper we attempt to answer the first two of 
those questions. 



II. RELAXATION OF THE LONG-RANGE 
FLUCTUATIONS 
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A. The setting 

In this section we develop analytical description of 
the time evolution of fluctuations due to rescattering of 
charged particles in the hadronic gas at the late stages of 
the heavy ion collision. [] The following physical picture 
underlines our approach. 

We consider distributions in rapidity space of a con- 
served net charge (e.g., electric or baryon). We start our 
description from the beginning of the hadronic phase. In 
each event this phase "inherits" a certain distribution in 
rapidity space of the charge from the primordial quark- 
gluon plasma phase. In the spirit of the statistical ap- 
proach, we do not consider a specific charge distribution, 
but an ensemble of distributions corresponding to a set 
of events. This ensemble of charge distributions is char- 
acterized by the mean and the magnitude of fluctuations. 
The initial magnitude of those fluctuations is the input 
of our calculation. 

The fluctuations of different length, or range, in ra- 
pidity space (i.e., different harmonics of the charge dis- 
tribution) relax on different time scales. The fact that 
the net charge is conserved is crucial here. Because the 
relaxation can only proceed via diffusion of the charge, 
the longer range fluctuations relax slower. The relax- 
ation time grows as a square of the range. Our goal is to 
provide a quantitative description of this process. 

In particular, we shall derive the equation which solves 
the following problem: Given a certain detector rapid- 
ity acceptance window, and given the initial magnitude 
of the charge fluctuations, to find the time evolution of 
that magnitude. From the fact that longer range charge 
fluctuations relax slower, it is evident that fluctuations 
of total charge in a wider rapidity window relax slower. 
Using typical numbers for the life-time of the hadronic 
phase and rescattering properties of the hadronic gas we 
shall determine how wide the rapidity window must be 
chosen in order to preserve, to a sufficient extent, the 
primordial (small) fluctuation magnitude. 



B. Diffusion in rapidity space 

We introduce the density of a conserved net charge Q 
per unit rapidity, n(y,r), as a function of time r in a 
given event. The objects we are really interested in are 
the fluctuations of the density, 



f{y,r) = n{y,r) - no, 



(1) 



* Our approach is different from that of |lo|-|l5| . We provide 
quantitative description of the time evolution of long-range 
fluctuations based on a stochastic diffusion equation, which 
we solve analytically. The value of the diffusion coefficient is 
our only numerical input and we determine it using realistic 
hadronic gas properties. 



In the Bjorken boost invariant scenario no is propor- 
tional to the total charge in the central region, no — 
Qtotai/ytotai, which is a constant within each event. 

We begin by considering evolution of / in a given event. 
The event - by-e vent fluctuations will be the subject of 



subsection II D . There we shall also address fluctuations 
of the total charge (Q, or no), which we ignore for now. 
Since / is a density of a conserved quantity it should obey 
a diffusion equation in rapidity space: 

d T f = ld 2 y f, (2) 

where 7 = 7(7-) is the diffusion coefficient. This diffusion 
equation (|2|) applies, provided the following conditions 
are satisfied: 

(i) We must assume that fluctuations are sufficiently 
small, to justify linear approximation in (^). This is rea- 
sonable if the number of particles carrying the charge is 
large, which is fulfilled in sufficiently central heavy ion 
collisions at the SPS and higher energies. 

(ii) We must assume separation of scales — the minimal 
interval of rapidity we can consider must be much larger 
than the mean rapidity change of a charged particle in 
a collision, 8y co \\. In other words, / must be averaged 
over sufficiently wide interval Ay. As we shall see later, 
the typical <5y co n for the electrical charge is of order 0.8, 
while for the baryon charge it is 0.2. Even in such a large 
acceptance detector as STAR at RHIC, the coverage does 
not exceed few units of rapidity. However, we can still 
consider a somewhat idealized limit of rapidity windows 
much wider than 5y co \\. This can certainly be reasonable 
for the diffusion of the baryon charge. Moreover, for the 
ideas we are evaluating to work, this limit is a necessary 
condition |l^,|TI| . Indeed, if the rapidity window is not 
large compared to Sy co \\, the fluctuations will relax by 
ballistic transport of the charge on the time scale of the 
mean free time Tf reo towards their thermodynamic values, 
too fast for those ideas to work. [[] 

It will not be essential for the analysis of this section 
how the coefficient 7(7-) is calculated. It may be (ide- 
ally) determined from first principles of QCD. In the late 
hadronic stage of the collision, which is our concern in 
this paper, one can also use a simple formula to estimate 
7. Since, on the microscopic level, 7 is due to particle 
collisions: [j] 

(3 ) 



2r fr , 



^In the same spirit, we shall not consider also the effects of 
finiteness of y ma x, the total collision rapidity range. These 
have been studied already, |I^,[l3|, and correspond to certain 
boundary conditions in equation (^). 

One can think of this as a random walk in rapidity space 
with Gaussian random steps of mean square length (5y co n) 2 
and compare the equation for the mean square distance from 
origin: (Ay) 2 = (r/rf rco )(5y co n) 2 with the Gaussian solution 
of the diffusion equation (0), which gives: (Ay) 2 = 2jt. 
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Since we shall be using Eq.(g) in Section [Fn] to estimate 
7 (or, rather, its time integral, see below), we shall discuss 
it here. It is clear that in the collisionless limit, Tf ree — > 
oo, the diffision in rapidity is absent, which is in trivial 
agreement with Eq.([|). It is instructive to consider also 
the limit of Tf ree ^ — the ideal hydrodynamic limit. In 
this limit diffusion is also frozen — particles have no time 
to propagate. Eq. (^) can still be applied, but one must 
realize that, effectively, 6y co n also vanishes in this limit. 
The reason is that successive rapidity shifts (each of order 
unity) are not independent but strongly anticorrelated. 
This happens because a colliding particle has no time to 
move out of the region of particles of its initial rapidity, 
and is more likely to scatter back in than further out. 
This is also related to the fact that space-time rapidity 
(1/2) ln[(t+z)/(t — z)] and momentum rapidity can differ 
significantly in this limit. 

In this paper we need to describe late not-so-dense 
hadronic stage of the collision with Tf rco only few times 
smaller than evolution time, and we shall apply eq.(||) 
without worrying about complications of the ideal hy- 
dro limit. One can also use cascade codes to estimate 
7, or just the diffusion length squared 2 J jdr, directly, 
with possible anticorrelations already included, bypass- 
ing evaluation of 5y co n- 



C. Langevin equation 

Equation (||) describes relaxation of the density n to 
its equilibrium value. It implies that we average over 
time scales longer compared to the characteristic time of 
fluctuation of /. If, however, we are interested in fluctu- 
ations of n around its equilibrium value, equation (^) is 
not enough. We must add a noise term: 



(f(y)f(y'))=x(y-y'), 



(6) 



(4) 



The characteristic autocorrelation time of the noise is of 
the same order as the time scale of the variation of /. 
This time is very short since the rapidity windows we 
consider contain many particles: / changes each time 
any of the particles collides. 

The rapidity correlation of the noise can be determined 
by requiring that equilibrium fluctuations of / are given 
by thermodynamics. Thus, they are Gaussian with the 
probability distribution 



J/] 



exp 



dydy' 



mm 

2x(y-y') 



(5) 



where x 1S the equilibrium density-density correlation 
function in rapidity, which we need not know explicitly 
at this pointj|, and we used boost invariance of Bjorken 
scenario. In other words: 



where the average is over time. Writing the functional 
Fokker-Plank equation for the probability density P[f] 
following from the Langevin equation (|J) we find that 
(||) is a stationary solution when the noise £ is Gaussian 
and obeys: 

T)Z(y', r0) - -2 7 a*x(y - y') 8{t - tO. (7) 



D. Fokker-Plank equation 

Before we proceed to solving the Langevin equation 
([|) , we should relate the time average in (|j) to the event- 
by-event average. They are the same by ergodicity. How- 
ever, one should also keep in mind that changes in initial 
conditions will change the equilibrium value n . Such 
initial state fluctuations, as discussed in the introduc- 
tion, can contribute significantly to the observed event- 
by-event fluctuations of the extensive quantities, such 
as the total charge Q. Such fluctuations are statisti- 
cally independent from the thermodynamic fluctuations 
we consider and can be calculated separately and added 
in quadratures, as is done in In contrast, for intensive 
quantities, such as, for example, the ratio of positively to 
negatively charged particles considered in R] , such initial 
state fluctuations cancel and need not be considered. 

Our goal now is to use equation (Q) to determine how 
quickly fluctuations approach their thermodynamic dis- 
tribution (||). In other words, we need to solve the func- 
tional Fokker-Plank equation with a given initial distri- 
bution P[f], and see how this functional evolves with 
time. 

The Fokker-Plank equation corresponding to (|^), (]?]) 
is given by 



d T P = dydy' 7 



My) 



- d 2 y f(y)P 
dlxiy - 2/0 



p 



(8) 



and is not very easy to study. Boost invariance helps, 
however. We shall do a Fourier transform with respect 
to the variable y. Let us denote harmonics of f(y) by 
fk- If we start with the factorized probability distribu- 
tion P = rifc-ffe [fk]> then each harmonic evolves inde- 
pendently following partial differential equation 



d 



f k Pk 



Xk-zrj-Pk 
dfk 



(9) 



It is easier to derive this equation directly from the 
Fourier transformed Langevin equation p*| 



^Relaxation time of small linear fluctuations we study will 
not depend on their absolute magnitude. 



*The condition of applicability of the diffusion approxima- 
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d T fk 



(10) 



E. Gaussian solution 

We can now pose the problem mathematically and 
solve it using equation (^). Let us imagine that at some 
initial time r = (in our analysis - the beginning of 
the hadronic phase) the probability distribution for each 
harmonic f k is given by a Gaussian with a mean square 
variance (fl) = cr 2 . What happens to this probabil- 
ity distribution with time? Substituting the Gaussian 
Ansatz 



P fe [/,] = (2^ (T)r i/2 exp J 



fl 



(11) 



one sees that equation (||) preserves the Gaussian shape 
of the probability distribution. The evolution of the dis- 
tribution can therefore be described by the evolution of 
crfc(T). The equation for <Jkij) following from (^) is given 
by 



Or 



with the solution: 



-2 7 fc 2 (a 2 - Xfe ), 



(12) 



4(r) =Xfc + (^(0)-X fc )exp|-2fc 2 y o jdrj. (13) 

For simplicity, we neglected the dependence of \k on T 
in the hadronic phase. Obviously, eq. ( p"2| ) can be solved 
for arbitrary Xfc( T )j if one is given. In this paper we 
apply eq. (12) to an idealized scenario, in which Xk 
changes abruptly at the QGP-to-hadron-gas transition, 
and see how the charge fluctuation magnitude relaxes to 
its hadron gas equilibrium value. Eq. (|13J) demonstrates 
explicitly that, due to the local conservation of charge, 
harmonics with small k relax very slowly. Using eq. (0), 
the exponent in (13) can be also written as k 2 (Aydis) 2 , 
where 



(A 2/diff ) 2 = 2 / 1 dr = 

Jo Jo 



= / {SycoiiY 



dr_ 

Tfrcc 



(14) 



is the mean distance on which a charged particle diffuses 
during the time r. 



tion, already discussed above, can be also stated as the condi- 
tion that we are describing sufficiently long-range, i.e., small 
k harmonics: k < l/Sy co n. For high harmonics the relaxation 
rate is l/rf rcc instead of -yk 2 . 



F. Relaxation in a rapidity window 

In experiment one is measuring the event-by-event fluc- 
tuation of charge, 



ndy, 



(15) 



in an interval of rapidity, [y\ , j/2] ■ This fluctuation, 
AQ = Q — (Q), can be written as a weighted sum of 
the fluctuations of the Fourier components /&: 



((AQ) 2 ) = 



dk 
2^ 



Ck? (fl) 



dk 
2^ 



44(r). (16) 



The weight is determined by the Fourier harmonics, c^, 
of the function equal to unity for y S [yi , 2/2] and zero 
otherwise, i.e., 



sin 2 (fcAy/2) 
(fc/2) 2 



(Ay = y 2 - yi ). (17) 



The main contribution is from harmonics with k ~ 1/Aj/ 
and smaller. Equations (|l3|) and ([l6]) determine the evo- 
lution of the magnitude of the fluctuation ((AQ) 2 ). 

For example, consider the case when the hadron gas 
"inherits" from the QGP phase practically zero fluctu- 
ations of the charge (or fluctuations much smaller than 
require d by hadron gas thermodynamics — an idealiza- 
tion of MP)), i.e., a 2 (0) « X k- Then, using eqs. ©, 
(|l4|), (|16|), and (|l7|), we can calculate the magnitude of 
the fluctuation of the charge Q in the rapidity window 
Ay after time r: 



«AQn = - 



dk sin 2 (fcAy/2) 
27 



(fc/2) 2 

Xk (l - exp {-fc 2 (Ay diff ) 
XoAyF(x), 

Ay d iff 
" (Ay/2)' 



! }) 



(18) 



The dependence on time r and rapidity window Ay en- 
ters through the ratio x of Audig given by (|TJ) and half 
width of the window Ay/2. [J The function F is given 
by: 



^)=l + ^(l-e-V^)-erf0) 



(19) 



It rises as x/^/tt at small x and saturates as 1 — 1/(x^/tt) 
at very large x (Fig.Q). xoAy is the equilibrium thermo- 
dynamic size of the fluctuation ((AQ) 2 ) (it scales linearly 



^We assumed, for simplicity, that Xk = const in the interval 
k < 1/ Ay, which is true when the density-density correlation 
length is much smaller than Ay — a good approximation for 
the large rapidity windows we consider: Ay ^> 8y co \\. 
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with the size Ay as it should). The reason that F satu- 
rates only as a power of time, and not exponentially, is 
the fact that characteristic relaxation time of the long- 
range harmonics diverges. 




0.5 1 1.5 



FIG. 1. The function describing onset and saturation of 
the equilibrium fluctuation magnitude in a rapidity window 
with time: (jT|), (|l|). 



III. REALISTIC RESCATTERING IN 
HADRONIC GAS 

A. Cascade estimates of the diffusion rates 

Our estimates for Az/diff in the hadronic gas phase can 
be simplified by noting that the dependence of 7 on r 
in (jlj) is mostly due to Tf reo . The rapidity change per 
collision, 5y co u, is approximately independent of r, since 
the decrease of the temperature is relatively weak in the 
whole hadronic stage, T = 160 — 110 MeV, and kine- 
matics of low energy scattering changes little. In other 
words, Aydiff is given by the random walk formula 



where 



dr 

3 

Tree 



(20) 



(21) 



is the total number of collisions per particle during the 
time r, from the beginning of the hadronic phase, r = 0, 
until freeze-out. 

At comparatively low relative energies of re-scattering 
in question, all relevant processes are due to certain res- 
onances, such as A for irN collisions, p for irir ones, etc. 
So, our first task is to determine what is the mean change 
of rapidity in such collisions. Taking pions and nucleons 
with thermal distribution, we have found probabilities of 
different resonances and have determined the following 
values for the mean rapidity shifts per collision 



The second (and much more involved) task is to eval- 
uate the average number of re-scatterings, for each type 
of particles. The average value can be evaluated from 
standard rates 



E 



dt crij Vij(t) n,j(t) 



(23) 



and even a very simple estimate shows that those num- 
bers are not small. The nucleons suffer especially multi- 
ple collisions, of the order of 10-20, since the A-induced 
cross section at its peak is as large as 200 mb, and there 
are many pions around. (This phenomenon is sometimes 
referred to as pion wind.) However, more accurate esti- 
mates depend on many details, like the geometrical loca- 
tion of the particle in question, expansion of the fireball 
etc. Those can only be obtained from either realistic hy- 
dro simulations or from hadronic cascades. 

We have chosen the latter alternative, using two pop- 
ular hadronic cascade codes. First, we have generated 
AuAu RHIC events (100+100 GeV) with RQMD jjj), 
and looked at the total rescattering number per particle 
located at midrapidity y rs 0. The results are shown in 
Fig.^, as a distribution over the number of rescatterings 
suffered by Kaons and Nucleons. One can see that for the 
nucleons the distribution has a very wide tail, reaching 
large values and with an average being of the order of 20. 
Kaons, on the other hand, suffer much smaller number 
of collisions, about 5 in average, and there is about 20% 
of kaons without rescattering at all. 

Unfortunately standard RQMD output does not allow 
to trace an individual pion, which appears and disap- 
pears, and we do not have the mean rescattering pion 
number directly from that simulation. We can only eval- 
uate the total number of n — N rescatterings 



N v — N N 
Jv coll — Jv coll 



(dN N /dy) 
(dN^/dy) 



(24) 



but those are far from being dominant. 

Fortunately, a similar code, UrQMD p7|, has been ex- 
tensively studied, and in the talk by M. Bleicher one finds 
a graph showing dN co n/ds for baryon-baryon, baryon- 
meson and meson-meson collisions, under the same con- 
ditions. Dividing by the number of mesons and integrat- 
ing over subcollision invariant energy, s, we obtained the 
total number of collisions from meson-meson scattering 
to be N™f on ss 8 per meson. Most of the mesons are pi- 
ons, of course, and this number fits well between N and 
K ones mentioned above. 



coll 



0.. 



N 
coll 



0.2. 



(22) 
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Number of rescatterings 

FIG. 2. Distribution over the number of rescatterings in 
central AuAu 100+100 AGeV collisions, according to RQMD 
code. The solid line is for Kaons, the dashed one for Nucleons, 
both taken in the mid-rapidity region (\y\ < 2). 

Adopting these numbers, we can combine it with ra- 
pidity change per scattering, and get our final estimates 
for the diffusion distance at RHIC: 



A^ iff « 2.2; 



0.9. 



(25) 



At SPS the pion multiplicity is about factor 2 smaller, 
and the rescattering number is reduced to N^ n w 15 and 



/V Tr 
JV coll 



5, so that at SPS 



Ay^ 



1. 



Ay, 



AT 
diff 



0.77, 



(26) 



Let us now compare it to some typical experimental con- 
ditions at SPS. For the accepted interval of rapidity of 
Ay = 2, the argument of F(x) plotted in Fig.[l] is di- 
rectly these values of Ay as, and one can read F from it: 
the result is 0.7 and 0.42 for pions and nucleons. 

In the previous section we described an idealized case 
c 2 (0) -C Xk- For our estimate of the effect we shall 
use more general expression, which gives the time depen- 
dence of ((AQ) 2 ) given its initial value at the beginning 
of the hadronic phase, r = 0, ((AQ) 2 ) = cr 2 (0)Ay [TJand 



the hadronic gas equilibrium value ((AQY 



XoAy. 



<(AQ) 2 



((AQ) 2 ) oq F(x) 

_ Aj/dig 
X ~ (Ay/2)' 



((Ag) 2 )o (1 



F(x)); 
(27) 



For example, starting from the fluctuations suppressed 
by a factor of about 3 i.e., ((AQ) 2 ) /((AQ) 2 ) eq w 

1/3, and using F » 0.70 for the fluctuations of the electric 



charge we find: ((AQ) 2 )/((AQ) 2 



'cq 



0.80. I.e., only a 



20% suppression survives until freeze-out after a 3-fold 
suppression at the QGP-hadron-gas transition. Similar 
estimate for the baryon number with F w 0.42, is more 



encouraging: ((AS) 2 )/ ((AS)' 



0.62. 



B. Charge fluctuations at SPS 

We have shown above, that the pion diffusion in ra- 
pidity in the hadronic phase is sufficiently strong. For a 
typical detector with the acceptance Ay = 2 it reduces 
the initial deviation from equilibrium value of event-by- 
event charge fluctuation to nearly its equilibrium value 
in the hadronic gas given by 



((AQ) 2 



0.7, 



ch 



(28) 



which is very far from the QGP value. The deviation 
from the Poissonian value of 1 is due to the 7r + , tt~ mul- 
tiplicity correlation induced by resonance decays. 

Furthermore, a specific centrality dependence of charge 
fluctuations in PbPb collisions follows from the resonance 
gas freezeout scenario. The magnitude of charge fluctua- 
tions (|2^) should increase towards most central collisions. 
This is opposite to what the QGP signal is supposed to 
do, if one expects it to show up in more central collisions. 

This opposite trend is due to the fact that more cen- 
tral collisions correspond to lower freeze-out temperature 
]l8[ : the larger multiplicity is, the later freeze-out oc- 
curs. Significant centrality dependence of the radial flow 
observed at SPS confirms this idea. It follows, that for 
most central collisions one expects less resonances, and 
correspondingly decreased deviation of charge fluctuation 
from the Poissonian value 1. 

A simple estimate of this effect can be made using the 
fact that the contribution of the resonances is controlled 
by the Boltzmann factor. Thus, for example, a decrease 
of temperature by AT w 20 MeV leads to a decrease 
of the contribution of the po resonance by a factor of 
2. Similar decrease should occur in the contribution of 
other resonances, such as u>. Therefore the variation with 
centrality of the freezeout temperature of order 20 MeV 
would correspond to centrality dependence of (|2|) of the 
order of 10%, rising towards more central events. 



"As with Xk> we neglect fc-dependence in at in the interval 
of k ~ 1/Ay. 
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C. Charge fluctuations at RHIC 

Finally, let us briefly discuss prospects of observation 
of QGP-like charge fluctuations at RHIC. STAR detec- 
tor has significantly larger rapidity acceptance, and al- 
though diffusion of pions is slightly stronger at RHIC, 
one can hope to see deviations from the resonance gas 
value. Indeed, for the rapidity window Ay = 4, equa- 
tions (pq), (H^l) predict that the primordial QGP factor 
1/3 suppression of electric charge fluctuation will survive 
as factor 



((AQ) ; 
<(AQ) 2 ) 



0.7 



(29) 



suppression. Another important prediction is that this 
suppression must strengthen, i.e., the ratio (29) decrease, 
as the width Ay is increased. In other words, opening 
wider rapidity "aperture" allows to see deeper back into 
the history of the collision. 

The baryon number fluctuations have slower diffusion, 
but rapid proton-neutron conversion together with vir- 
tual invisibility of neutrons basically undermine the very 
idea: the relaxation slow-down can only occur if the 
current under consideration is a conserved one. Similar 
problem arises if one considers strangeness fluctuations. 



IV. TT+/-K- FLUCTUATIONS NEAR THE 
CRITICAL POINT E 



In this section we discuss a new signature of the critical 



end-point E 119 
discussed in [H 



based on the tt + /tt~ fluctuations. As 
the main feature of thermodynamics 
near the critical point is the presence of long-wavelength 
fluctuations of the magnitude of the chiral condensate^ 
or, the cr-field. One of the signatures proposed in [^9|j4| 
was based on the fact that such light cr-quanta cannot 
decay at freeze-out, if it occurs near point E, and their 
large population survives until the time after freeze-out 
when their rising mass exceeds the 2ir threshold. The 
signature proposed in |l!|f| used the fact that produced 
pions had a necessarily non-thermal spectrum with low 
mean px ■ 

The contribution of a decays to the pion spectrum is, 
to some extent, similar to contributions of other, stan- 
dard, resonances in the resonance gas model, such as p°. 
In particular, the produced 7r + 7r~ pairs introduce posi- 
tive correlation between the numbers of positive and neg- 
ative pions, therefore reducing the fluctuation of n + /?t~ 



ratio. The most significant difference is that the contri- 
bution of a decays is only present near the critical point. 

In order to estimate this contribution, we assume that 
the thermal population of cr's is approximately half of the 
so-called "direct" charged pions (since the mass of the 
sigma is comparable to m n ), as is done in Q|. Thus the 
decaying cr's produce the number of charged pions of or- 
der 2x2/3x1/2 = 2/3 of the number of "direct" charged 
pions. The latter is about 1/3 of the total number, the 
rest are produced by resonance decays. Therefore, the 
contribution of sigmas to the charge pion multiplicity is 
about 2/3 X 1/3 w 20%. Thus, the contribution of cr's 
can be estimated to result in a reduction of fluctuations 
by an amount of order 20%. * 



Even though the actual contribution will also depend 
on the detector acceptance of the soft part of the pion 
spectrum, it is reasonable to expect a noticeable suppres- 
sion, comparable to the suppression of about 30% due to 
standard resonances (p°, uj, etc.) in (^8|). 

What is important is that such an additional suppres- 
sion of the TT + /n~ fluctuations can only appear if the 
freeze-out occurs near the critical point. As a function 
of the collision energy, this effect will be seen as a dip 
in the magnitude of the n + /ir~ fluctuation, thus provid- 
ing a signature for the discovery of E. Away from the 
critical point the fluctuations will be compatible with the 
ordinary resonance gas result calculated in . 

The NA49 collaboration have now acquired data from 
collisions at 40, 80 and 160 AGeV of initial energy at 
SPS. We are looking forward to the analysis of this data. 
If it shows a monotonic dependence of ir + /ir~ ratio fluc- 
tuation on collision energy as well as on centrality, consis- 
tent with the resonance gas prediction, then one should 
be able to rule out the presence of the critical point in 
regions of the QCD phase diagram the location of which 
can be determined by the analysis similar to J2(| ■ Colli- 
sions at other energies, in particular those of RHIC, will 
be needed to scan other regions of the phase diagram. If 
a signal described here is observed, other signatures dis- 
cussed in Q must also be seen to confirm the discovery 
of the critical point. 
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§§ As emphasized in |hJ-[L3|, Ay should be sufficiently smaller 
than the total collision rapidity range j/ max . Otherwise the 
fluctuations will be trivially constrained by the total charge 
conservation. In a sense, we would be looking beyond the 
QGP stage of the collision. 



*** This can be compared to the reduction due to p° res- 
onance, which is about 9% @. Another 20% m eq.(|Z8|) is 
contributed by u) and other heavier resonances. 
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